The Interplay of Charge and Spin in Quantum Dots: The Ising Case 
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The physics of quantum dots is succinctly depicted by the Universal Hamiltonian, where only zero 
mode interactions are included. In the case where the latter involve charging and isotropic spin- 
exchange terms, this would lead to a non-Abelian action. Here we address an Ising spin-exchange 
interaction, which leads to an Abelian action. The analysis of this simplified yet non-trivial model 
shed some light on a more general case of charge and spin entanglement. We present a calculation of 
the tunneling density of states and of the dynamic magnetic susceptibility. Our results are amenable 
to experimental study and may allow for an experimental determination of the exchange interaction 
strength. 
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I. INTRODUCTION 

Significant progress in the study of the physics of quan- 
tum dots (QDs) has been achieved following the intro- 
duction of the Universal HamiltoniarJn^ (UH). The latter 
facilitated the simplification of intricate electron-electron 
interactions within a QD in a controlled way. Within that 
scheme interactions are represented as the sum of three 
spatially independent terms: charging, spin- exchange, 
and Cooper channel. Notably, even the inclusion of the 
first two terms turned out to be non-trivial: the resulting 
action is non-Abelian^. 

To understand the complexity of such a problem one 
can refer to the case of eharging-only interaction. As 
was suggested by Kamenev and Geieo^, one can take 
the following steps in solving that problem: start from 
a fermionic action which includes an interaction term 
quadric in the (fermionic Grassman) variables, perform a 
Hubbard-Stratonovich transformation by introducing an 
auxiliary bosonic field, then perform a gauge transfor- 
mation over the Grassman variables, and finally integrate 
them out. The resulting, purely bosonic, action is simple. 
In an imaginary time (Matsubara) picture the action is 
quadratic in the bosonic components, which renders this 
action easily solvable. The trick of gauge-integrating over 
Grassman variables does not work for the non-Abelian 
case^ so that an alternative approach is needed. 

Attempts to account for charge and spin interactions 
in QD have been reported earlier. Alhassid and Rupp^ 
have found an exact solution for the partition function 
(and susceptibility); elements of their analysis were then 
incorporated in a master equation analysis of transport 
through the QD. More recently an exact solution of the 
isotropic spin interaction model has been presented^. For 
the latter model some quantities turn out to be particu- 
larly simple (e.g. the finite frequency spin susceptibility 
vanishes; evidently there is no difference between longi- 
tudinal and transverse spin susceptibility). This means 
that the analysis of a model with anisotropy in the spin 
interaction is called for. A perturbation expansion in spin 



anisotropy has been reported earlier—, but it still remains 
desirable to consider an anisotropic model which can be 
analyzed exactly. By considering such a model one would 
be able to understand the entanglement between charge 
and spin degrees of freedom, and also see in detail how 
a non- vanishing, complex spin susceptibility arises. This 
is the focal point of the present analysis. 

In bulk systems the exchange interaction competes 
with the kinetic energy leading to Stoner Instability 
(SI)i^ In finite size systems mesoscopie Stoner unstable 
regime may be a precursor of bulk thermodynamic SI. 
We consider here an Ising spin interaction. Such a model 
is Abelian, and complications due to non-commutativity 
of different terms in the action do not arise here. Also 
such a model does not exhibit a mesoscopie Stoner un- 
stable regime^. This means that at zero temperature, as 
the dimensionless parameter J/ A (J being the exchange 
interaction strength and A is the mean level spacing) the 
system abruptly switches from a paramagnetic to a (ther- 
modynamic Stoner unstable) ferromagnetic phase. We 
stress that notwithstanding the simplicity of the model 
considered, spin-charge entanglement is present here, and 
non-trivial transverse a.e. susceptibility does arise. Some 
of our conclusions can in principle tested in QDs made of 
materials close to the thermodynamic Stoner Instability, 
e.g., Co impurities in Pd or Pt host, Fe or Mn dissolved in 
various transition metal alloys, Ni impurities in Pd host, 
and Co in Fe grains, as well as new nearly ferromagnetic 
rare earth materials A"— 

The outline of this paper is as follows. In Section II 
we introduce our model Hamiltonian and the subsequent 
imaginary time action. In Section III we employ the 
technique of zero-dimensional functional bosonization^ 
which eventually allows us to express the single-particle 
Green's Function as a product of the non-interacting 
Green's Function and a term which depends on two 
bosonic fields. We then show how to reduce the problem 
to that of classical stochastic equations for the bosonic 
fields. In section IV we express the grand-canonical par- 
tition function in terms of canonical ones, leading to both 



a mathematical and physical sim.plification of the calcu- 
lation. In Section V we calculate the tunneling density of 
states and in Section VI longitudinal and transverse spin 
susceptibilities. Section VII presents a summary of the 
main results with some perspectives. We include some 
more technical calculations in three Appendices. 



II. HAMILTONIAN AND EFFECTIVE ACTION 

We consider a normal-metal QD in the metallic regime, 
where the Thouless energy Etii and the mean level spac- 
ing A satisfy g = Et^/A » 1 (g is the dimensionless 
conductance) and a temperature T 3> A. It is the regime 
where a description in terms of UH is viable. 

We restrict ourselves to a simplified version of the UH 
where the interaction in the Cooper channel is set to 
zero and the spin-exchange term is chosen to be a fully 
anisotropic Ising-like term, —JS^, with a ferromagnetic 
exchange coupling, J > 0, Sz is the total spin of the dot 
in the z direction. This form of interaction is sufficient 
to bring about the Stoner instability phenomenon and 
other spin-related effects, whilst avoiding calculational 
complexities inherent to a fully spin-symmetric model. 
Possible physical sources for such an anisotropy may 
include geometrical and/or molecular anisotropy, mag- 
netic impurities in the system, or even the application of 
anisotropic mechanical pressure. 

The complete form of the reduced UH is thus 

a.,(T a.(T 

- -^[^aLa<^aa'aa^a'\ ■ (1) 
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Here {Sa} is a set of electronic levels in the dot, and A^o 
in the charging term represents a positive background 
charge controlled via an external gate. We assume that 
the QD is either isolated or weakly-coupled to the leads 
and in the Coulomb blockade regime. On the other hand, 
we will be considering the spin-disordered regime below 
the Stoner instability. So the parameters of the Hamil- 
tonian ([T]) obey 



J < A < T < £;c 



(2) 



where T = /3~^ is the temperature. 

The Euclidean action corresponding to the Hamilto- 
nian ([T]) is given by 

5[¥, *] = 51 / dr j^c<(a. + e. - M)*a 



Sc[5I*"*"-^o]'-^[^^a(T"^a]'|, 



(3) 



where we use spinor notations ^^ = (■0-1'q(t), V'^q ('''))■ 
We introduce two auxiliary bosonic fields, <^'^(t) and 



f^ir), to decouple the Coulomb and exchange terms with 
the help of a standard Hubbard-Stratonovich (HS) trans- 
formation. This results in the following action: 



5 = 5"= 



(4) 



where 
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5"== / dr 
Jo 

5-^ = / dr-^--^ 
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(5) 



S'^'^ = / dr V *„ [9, + e„ - /i + i^- + a'ip'] *„. 

Here a^ is a Pauli matrix, the bosonic fields are periodic 
and the fermionic fields are antiperiodic in r with period 
/3. This action is the starting point for all the subsequent 
calculations. We will use the functional bosonization ap- 
proach as developed in^— : first we gauge out the mixed 
fermionic-bosonic terms in the action ([5|) and then inte- 
grate over the fermionic field thus arriving at a purely 
bosonic action. After that, instead of dealing with this 
action directly we will use a stochastic bosonization as 
described in the following section. 



III. FROM FUNCTIONAL TO STOCHASTIC 
BOSONIZATION 

In order to gauge out the mixed fermionic-bosonic 
terms in the action ^ , we introduce a generalized gauge 

transformation, ^^ = T^^'^a , ^a = 'J'q.'J with 



j,i0'^(r)+r(r) 

e"'"^^)-*'^^) 



"Gauging out" implies the following identity 

*„ [dr + i(^^(r) + a'ip^T)] *„ - %^ [dr+A] *„ , (6) 

where A is some constant matrix. In order to fulfill (|6]) 
we require the gauge matrix T to obey 

[dr + iip'{T) + a'ip''{T)]7^7A. (7) 

Since the bosonic fields are real, this equation sepa- 
rates into real and imaginary parts, corresponding to the 
exchange and charge channels. Using the substitution 
A = A^a^ -I- iyi'^ for the constant matrix A in the matrix 
gauge equation ([7|), we have 



9°(r)==yi''-(^"(r) 



(8) 



where a stands either for charge, c, or for spin, s. 

To determine the constants ^l*^ and A'^ we note that 
the antipcriodicity of the fermionic fields requires that 
T(/3) = T(0). This in turn implies 61" (/3) = d\^) 4- 27rin" 



and 9" (13) = 9''{0) + 2TrN with integer n" and N. Now we 
single out zero-Matsubara-frequency components of the 
bosonic fields (P°'{t): 



(^-(T)=y.» + ^«(r) 
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dTifi^ir). (9) 



Integrating Eqs. ([5]) over r from to /3 results in A'^ = 
ip^ + {2n/l3)N and ^l'' = (^^ + (27ri/;3)n" so that the gauge 
equations ^ reduce to the following form: 
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(10a) 
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After the gauge transformation the mixed action in 
Eq. ([5]) is reduced to the following quadratic fermionic 
action in terms of the transformed fields: 
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Sf= / dT^vl/„(r) 



dr +ea~~ fJ.0 
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(11) 



The zeroth components of the bosonic fields (|D5I) enters 
Eq. ()lip via the spin-dependent effective chemical poten- 
tial fia given by 

~ 2TTi 

Ha = l-i ~ i^o - (Tfl - ^ {N + an') , (12) 

where ct = ±1 for spin up/down respectively. 

The gauge equations (jlOp become important for corre- 
lation functions which are not gauge invariant but depend 
on phase terms which are functions of the gauge parame- 
ters 6'^ and 6' (e.g. the Green's function calculated in the 
following Section and Appendix |X]). These parameters 
are functionals of the bosonic fields (p'^{t) and ip^ir) re- 
spectively. Thus, in order to calculate these phase terms, 
one should solve the gauge equations and then carry out 
the integration over the bosonic fields4iiii. 

Here however, we consider an alternative method, 
which bypasses the need to carry out the functional in- 
tegrals over ip'^''{T). Even though in our case these inte- 
grations pose no great difficulty, the method we consider 
has general applicability and could be used in cases where 
such integrations are impossible to perform analytically. 

Our approach is to view the gauge equations ([TU]) as 
classical Langevin equations governing the stochastic dy- 
namics of 9'^ and 6', with the bosonic fields playing the 
role of noise. The distribution of the noise is determined 
by the bosonic actions S^ and S', Eq. ([S])- 

The Langevin equations can be mapped, via the stan- 
dard tools of classical stochastic analysis^, to Fokker- 
Planck (FP) equations from which the time dependent 
distribution functions for 9'^ and 9' can be determined. 
As an example, the form of the FP equation derived from 
Eq. ([TUa)) is 



dT'' /27r,^ .A aj"= 



E,- 



Q2yc 



(13) 



where 3"^ is the distribution function for the gauge pa- 
rameter 9'^ and C is a constant (details regarding the tran- 
sition from Langevin to FP equations and their solution 
are given in Appendix [C|) . Equation (jl3|) is a standard 
diffusion equation with a drift term, the solution of which 
(with an appropriate initial condition) is simply a decay- 
ing Gaussian, explicitly given by Eq. (jCSp . 

This distribution, and a similar one for 9s, can now 
be used to calculate the averaging of any phase terms 
involving the gauge parameters in the calculation of non 
gauge-invariant correlation functions. Thus we can, in 
effect, replace a functional integration with an integration 
over a finite number of parameters. This is an alternative 
method by which to integrate out the finite frequency 
components of the bosonic fields pf{T) and pf{T). 



IV. SINGLE PARTICLE GREEN'S FUNCTION: 
EFFECTIVE CHARGE QUANTIZATION 

We begin with calculating the temperature Green's 
function (GF) in the grand canonical ensemble and will 
show that in the Coulomb blockade regime it reduces 
naturally to one in the canonical ensemble. Our starting 
expression is: 



S<t(t, ^J) =^gQ,CT(r, /x) 



(14) 



g„., ^ _l^y D[vI/„vl/Je-^[^°*°]vI/„,,(r,)*a..(T/), 

where 9q,ct is an auxiliary GF corresponding to a level 
£q,, S [^Q^a] is the a-term in the Euclidean action ([3]) 
and T = Tf — Ti. 

After the HS transformation and gauge transform ([B|), 
the Gaussian integration over the quadratic fermionic 
action (fTTj) is straightforward. The resulting GF of 
non-interacting electrons corresponding to this action, 
Sa.ai'^^V-cr), depends - via Ea. ([T^ - only on the zero- 
frequency component of the bosonic fields (pg. This al- 
lows us to subdivide the remaining functional integration 
with the bosonic part of the action ([5]) into that over the 
zero-frequency, ipg, and finite frequency, ^^, components, 
which results in the following expression: 



Sa.a=n^(T)ff(r) 



((^°(M)S°,.(r,Ma)))o 
P"(m)))o 



(15) 



Here H''(r) are the phase correlation functions result- 
ing from the functional averaging of the charge or spin 
phase factors over the finite-frequency components of 
the appropriate fields, ((. . .))q stand for the functional 
integrals over the zeroth-component fields ipg and (Pq. 
All these functional integrals are defined in Eq. (|A1[) - 
(|X3l) in Appendix El Then Z°(/x) = Z'^{Jl^)Zl{ln) and 

Sq ct('''' A'cr) i^ 11^^ grand canonical partition functioni^ of 
non-interacting electrons with the spin-dependent chem- 
ical potential pt^, defined by Eq. ([T^ . 



The charging effects can be fully accounted for by in- 
troducing winding numbers in the integration over ipQ: 



where 



V'o 






27r 



(16) 



where — tt < ^q < tt and an integer tti is a winding num- 
ber. In the original work of Gefen and Kamene\i^ these 
were not considered, leading to an incorrect final result. 
They were first introduced in the context of the charg- 
ing interaction on small metallic grains by Efetov and 
Tschersicbi^ within a Matsubara framework, and were fi- 
nally correctly implemented by Sedlmayr, Yurkevich and 
Lerner— within a Keldysh-technique framework. The 
introduction of the winding numbers (J16p allows us to 
replace integration over Lp^ with summation over all in- 
tegers 771 and integration over Lp^. The sum over m is 
performed using the Poisson formula, which results in a 
new summation of the form 



^e-''^^(^-^°)'xJ(7V). 



N 

The Poisson resummation transforms summation over m 
into summation over the conjugate variable, N . In our 
case ipQ represents a phase, whose conjugate is evidently 
the particle number N . While the sum over the parame- 
ter m had many contributions (since {l3Ec)~^ ^ l)j the 
sum over N contains, under the conditions ([2]), only two 
terms N = A'q ± i, near the Coulomb peak {Nq is half an 
integer) and one term in the Coulomb valleys (i.e. every- 
where outside of the region of width T near the peak): 
the contribution of all the other terms is exponentially 
suppressed. This is a manifestation of charge quantiza- 
tion in QDs. 

In this way we perform the integration in Eq. (|15p to 
find (see Appendix [A]) : 

9a,^(T,M) = 2^e '' ^>> "^^Mat, (17) 

^W N 

CO 7T 

J.-/d^^e-^/^eK-+t)--ZO(^)gO^^(,,^,), 

— OO — TT 

(18) 

where the reduced phase correlation functions 11° are 
defined in Eq. (jA4p . The effective charge quantization 
in Eq. (|17p makes it natural to change over from grand 
canonical to canonical quantities for a given N, followed 
by a weighted summation over N, where required. Let us 
stress that the canonical quantities are auxiliary and we 
calculate in this way the grand canonical GF of Eq. ([T7| . 
Expressing 3n via canonical quantities leads to an ex- 
tra summation since 2, = ^^jC^^'^Zn, etc. This calcula- 
tion is detailed in Appendix [BJ The resulting full single 
particle GF in imaginary time (following summation over 
all single particle energy states) is given by 



S(r,/i) 



A sin(7r|T|T) F{0)' 



(19) 



N 



N M=-N 

SN = N-No-^^ En.m = 2EcSN-^. 
zhic z 

The double summation above arises from replacing the 
grand canonical partition function in terms of the sum 
over canonical ones, Z(/i) = X^n c^^^Zn- The summation 
parameters are the electron number, N , and the total 
spin of the dot (in the units of h/2), M . Naturally, the 
GF is spin independent: we are considering the regime of 
parameters, Eq. Q, below the Stoner instability where 
there is no symmetry breaking to distinguish opposite 
spin polarizations. Note that this result is valid in the 
regime ([2]), provided that 



iVA >r. 



(7V-|M|)A>T^ 



(21) 



i.e. when the QD contains many electrons and is not very 
close to the Stoner instability. Moreover, under these 
conditions the sum over M in Eqs. ([19]) and 120]) can 
be replaced by an integral from — cxd to +oo and the ex- 
ponent of ^gi^_ j\ resulting from this integration can be 
totally neglected. With the same accuracy, we should ne- 
glect the exchange energy J in the exponent in Eq. ()19p . 
Thus we find 



9ir,^l) 



1 



^Vc- 



/3Asin(^) Z^ 



^E^[p{5Nf~2TSN] /221 



N 



SO that under conditions © and (PT|) - not surprisingly - 
the one-particle GF is independent of the exchange part 
of the universal Hamiltonian ([1]). Such a dependence 
would emerge only very close to the Stoner instability, 
when |A — J|/ J <^ 1 but this parametric region is beyond 
the scope of the presented technique. 



V. TUNNELING DENSITY OF STATES 

The tunneling density of states (TDoS), i^(e), can 
be directly related to the conductance of the QD in 
the limit of weak coupling to the leads and is thus a 
quantity of great importance. The TDoS is given by 
iy{e) = -ilmg^(e), where the retarded GF, G^(e), is 
a Fourier transform of the GF in real time, G{t,^), ob- 
tained from Eq. ((22|) by the straightforward analytical 
continuation from the upper half-plane. Since g(r, /x) is 
independent of the exchange energy under the conditions 
da]) and Eq. ([21]), so is the TDoS.J^ 

For tutorial purposes, we use the results of Appen- 
dices E] and |B] to derive a more general expression for 
i/(e:), valid for any relation between the parameters in 
Eqs. ^ and (PT|) and show how it goes over to the known 
expressioni^ under conditions ^ and (PT|) . 



Using the GF in the e-rcprcscntation, Eq. (|B7|) . and 
performing the summation over all the levels as described 
at the end of Appendix |B] wc find 






1 ^ 

4y^ y^ g-/3Bc(N-JVo)'-i/3(A-J)A/=^ 

^ AT 71 T 7\r 



N M=-N 

1 - n{e ~ p.- S,n,m) +n{e 
where we have defined 



^jv-i,Af-i) , (23) 



^N,M = 2Ec{N ~ No + ^) - \J{M + i) , (24) 



and 



/2 = iA(7V + M) , 



No^No + ^, (25) 



while vq = 2/A is the TDoS in the absence of interac- 
tions, ji{e) = [l + c'^']~\ 

Equation ([25)) is the general expression for the TDoS 
for any combination of parameters for a many-electron 
dot. When the inequalities ([2]) and ([21]) are satisfied, 
we can easily sum over M as described at the end of 
the previous section and then limit the summation over 
N to the two terms for which the value of |A^ — 3\fo| is 
minimal (although deep in the Coulomb valley only one 
term is actually contributing). The resulting TDoS is 
independent of J (or, more precisely, tiny J-dependent 
corrections are beyond the accuracy of current calcula- 
tions and thus omitted) and coincides with that obtained 
in Ref. [ii 



iy{e) U{e-U) + c'f"-^''-'^^U{e~^N+i) 



vo 



1 -|- Q-P(iN-tA 



(26) 



where U{e — S,N) = 'n[e — ^N-i—li) + ^ — n{e — ^N — p), and 
S,N is obtained from £,n,m by putting J = in Eq. ([24)) . 
We illustrate the dependence of v on energy for integer, 
half-integer and intermediate values of iVo in Fig. [U for 
a specific choice of parameter values T and A. Its de- 
pendence on temperature at the bottom of a Coulomb 
blockade valley is depicted in Fig. O It is important to 
note that the TDoS obtained in the Coulomb valleys is 
not physical since we neglect co-tunneling contributions; 
however, the T-dependcncc near the peak will be ob- 
tained as a linear combinations of those shown in Fig. [21 
Note that for any given set of parameters the center of 

the TDoS curve is at Eq = ^ATVo - 2^c (n - Nq) and 

thus a function of A^o, Eq. (|25p . This moving from one 
Coulomb valley to the next, the TDoS curve is shifted by 
A/2 due to adding an extra electron to the dot, which 
raises the effective chemical potential and thus shift the 
TDoS curve. That is the reason for the 'half-gap' in 
TDoS at the degeneracy point. 



VI. MAGNETIC SUSCEPTIBILITY 

We now turn to calculating the longitudinal and trans- 
verse magnetic susceptibilities of the system. 
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FIG. 1. TDoS (in units of i^o) as a function of e = ejE^ for 
T = 0.2£;c and A/T = 0.1 in (a) a CB valley (iVo = 100), 
(b) an intermediate region (A^o = 100.35), (c) a CB peak 
(iVo = 100.5). 



It is clear that only the static component of the longi- 
tudinal susceptibility is non-zero due to the lack of spin 
flip processes in the Ising model.— A direct calculation 
of the correlation function (S'z(r)S'2(0)} shows this to be 
r-independent, as expected. The static susceptibility is 
given by 



1 d^ 

Xzz = -^ lim -—^\nZ{h), 



(27) 



where Z{h) is the partition function of the system calcu- 
lated in the presence of the following source term in the 




FIG. 2. Dependence of the TDoS (in units of z^o) on the 
temperature (measured in Ec) at the bottom of a CB valley 
(iVo = 70) for A/Ec = 0.02. 



action: 



^ "'0 



(28) 



The calculation is straightforward, leading to the result 
L{h) ~ K cxp ■^ 



\4/3(A-J)/' 

with K being some irrelevant constant. Plugging this into 
the definition ([27l) yields the well known expression 



Xzz{(^ = 0) = 



1 1 



2 A -J 



(29) 



As expected, the static susceptibility is independent of 
the number of particles on the dot, external gate voltage, 
charging effects, etc. 

We now turn to a calculation of the transverse mag- 
netic susceptibility. This quantity is inherently different 
from the longitudinal one since it is dynamic: the model 
allows for transitions between different transverse spin 
polarization states. 

We define the dynamic transverse susceptibility in 
imaginary time as 



X+-(r) = (a+(0)a-(r)). 



(30) 



where a+ = J^a *«T*"i and a' = J^a *ai*at- Thus 
we need to calculate the functional average of 

a,l3 

with the action given by Eq. ([3]). The procedure closely 
follows to that of the calculation of the GF described in 
Section IIVI The final outcome of this calculation is 

N 
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(31) 
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FIG. 3. -ilmx^ as a function of frequency lj (in units of 
A) for (a) A ^ 0.1 and f = 0.1,0.2,0.3 for the left, cen- 
ter and rightmost curves respectively and (b) -^ = 0.1 and 
A — 0.05,0.1,0.2 for the top, center and bottom curves re- 
spectively 



where ficr = N„A and N^^ is the total number of electrons 
with the spin projection a ="[, J,. 

Fourier-transforming the result of Eq. ([3T|) to Matsub- 
ara frequencies and then performing a simple analytic 
continuation, we find the imaginary part of the physical 
response function x^ (lu): 



Imx+-(.) = V^/^(A-J)^f[(A+.)-(A-.)^] 






sinh 


' I3iu' 
2 




sinh 


f3A 
2 


(1^ 


-7). 



(32) 



This function is depicted in Fig. [3] The most salient 
features are a linear dependence at the origin and the 
existence of a peak at a certain wq . Both the slope at the 
origin and the value of wq can be used used to charac- 
terize an experimentally obtained curve of the transverse 
magnetic susceptibility as a function of frequency. We 
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FIG. 4. Fit of numerically acquired data for FWHM to func- 
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aujo, yielding a = 1.59. R for this fit is 



find the slope at oj ~> as 
1 



-Imx+-(.^0)«^^^, 



(33) 



where the approximation was made consistent with the 
inequality Q. Under the same condition, the peak fre- 
quency is given by 



Wo 



2J2 



/3(A-J) 



(34) 



Yet another parameter of interest is the full-width-at- 
half- maximum (FWHM). Numerical analysis shows that 
it is proportional to the resonance frequency: FWHM k, 
1.59wo. This result was derived by numerically obtaining 
the FWHM for various values of wq and fitting the results 
to a linear curve, as shown in Fig. U] 

The imaginary part of the susceptibility represents the 
systems capacity to absorb and dissipate magnetic energy 
at a nonzero frequency. For the static susceptibility only 
a real part is finite. A simple calculation leads to 



Rcx^ 



0) 



cf(A+"') 
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(35) 



Note that in the limit ,/ = we recover the well known 
identity x^ — '^Xzz for the static susceptibilities. The 
real part of x^ at finite frequencies can be found either 
directly or via the Kramers-Kronig relations but we do 
not present the result here as it has little physical rele- 
vance. 

As in the case of the longitudinal magnetic susceptibil- 
ity, it is clear that the transverse susceptibility is not af- 
fected by the charging interaction in the dot. Once again 
we see that under conditions ([2]) and (piT) the charge and 
spin degrees of freedom are effectively decoupled. 



VII. SUMMARY 

The main results of this work fall into three basic cat- 
egories. These are the single particle GF, the TDoS, and 
the magnetic susceptibilities. The results for all three 
classes of correlation functions were obtained by means of 
the functional bosonization approach combined with the 
solution of classical stochastic equations for the bosonic 
fields. We considered the Ising version of the Universal 
Hamiltonian for description of the interplay between the 
spin and charge degrees of freedom in zero-dimensional 
systems. Such model is Abelian and therefore does not 
include the physics of non-commutative variables. It also 
does not exhibit the mesoscopic Stoner instability regime. 
Nevertheless, the spin-charge entanglement is present be- 
ing manifested in e.g. non-trivial AC spin susceptibil- 
ity. The model, being a simplified version of the quan- 
tum Universal Hamiltonian model gives qualitatively cor- 
rect description of the thermodynamics and transport 
through nanostructures in the vicinity of thermodynamic 
Stoner Instability point. The Stochastic Bosonization ap- 
pears to be very powerful tool for a treatment of Abelian 
gauge theories and a promising method for solving non- 
Abelian models corresponding to isotropic/anysotropic 
quantum limits of the Universal Hamiltonian. The the- 
ory of thermodynamic Stoner Instability and its influence 
on the transport through single electron transistor can be 
tested experimentally in quantum dot devices and gran- 
ular systems^"—. 

We summarize below the central results and key ob- 
servations reported in the paper. 

• Canonical variables and charge quantization. In 
our calculation of the GF, the tools we used and 
the choices made not only allowed us to carry out 
a non-perturbative calculation, but also had phys- 
ical significance. The use of functional bosoniza- 
tion and generalized gauge transformations and 
the implementation of winding numbers, as well as 
the transformation to conjugated variables via the 
Poisson rc-summation, led us to employ canonical 
quantities. The latter is a consequence of strong 
charging interaction. 

• Regimes of validity. The transition to canonical 
quantities, namely the introduction of the canon- 
ical partition function, also led to further insight 
with regard to the various physical regimes the 
system may be found in. Our calculation of the 
canonical partition function itself (and the associ- 
ated quantity 2,jv(^q)) imposed limitations on the 
physical parameters involved. We found that the 
system must be large enough (meaning a large num- 
ber of electrons), and far below the Stoner insta- 
bility point. We had to self consistently assume 
that the fiuctuations in the systems magnetization 
were much smaller than the system size. This cor- 
responds to a requirement that the system be far 



from a phase transition point, which in our case is 
the SI point. 



Appendix A: Grand Canonical Single Particle 
Green's Function 



• Spin- charge entanglement. Introduction of the 
canonical partition functions led directly to a sum- 
mation over all possible values of the magnetiza- 
tion. These are of course limited to \M\ < N. 
Since the number of particles itself is controlled by 
the charging interaction when in the CB regime, 
and the fluctuations of the magnetization are in- 
fluenced by the exchange interaction, this can be 
seen as a form of coupling between the charge and 
spin degrees of freedom. The coupling between the 
two interaction channels becomes important as the 
magnitude of magnetization fluctuations increases, 
i.e. as one approaches the SI point. Only then do 
values of M which approach the system size be- 
come accessible and, consequently, of physical im- 
portance. Far below the SI point, the spin-charge 
coupling is very weak, and effects of interplay are 
minimal. Our calculation of the TDoS showed the 
exchange interaction to have an extremely negli- 
gible effect. The magnetic susceptibilities in turn 
showed no dependence on the charging interaction. 



• Determining J and A. The calculation of the trans- 
verse magnetic susceptibility is, to our knowledge, 
a new result, and perhaps the most important in 
this work. As we have discussed previously, the 
importance of this result is that it provides an ex- 
perimental method to determine the values of the 
parameters J and A. Our result is a direct predic- 
tion of the absorption spectrum of the system, and 
as such should be amenable to experimental mea- 
surement. The various curve characteristics which 
we derived, including the slope at w ^' 0, the lo- 
cation of the resonance frequency and the FWHM, 
should in principal, through their dependence on J 
and A, allow these values to be ascertained from 
such a measurement. 
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In this appendix we present a detailed non- 
perturbative calculation of the single particle Green's 
Function (GF) for our model system ([T|). The GF it- 
self was used in order to derive the tunneling density of 
states (TDoS), but its calculation also serves to show the 
methodology used in calculating the various other quan- 
tities considered in this work. 

As discussed in section II, a HS transformation is 
applied, reducing the action to the form presented in 
Eqs. (HI and (O. Carrying out the Gaussian integration 
over the fermionic fields after the gauge transformation 
([TU]). we obtain the GF as follows: 

g„^,(r,/.)=ff(T)ff(r)x 
Jd^le-^<M,] /d^ge-^oK] [Z"(/I)go^^(r,/I,)] 

— OO — oo 

7 difi^o e-^ot"^"] 7 difil e-^°['^°] ZO(/I) 

— OO — OO 

(Al) 

Here Z°(/I) = Z'^{Jl^)Zl{]li) and Sa,a('r, A^<t) are the 
grand canonical partition functioiii^ and GF of non- 
interacting electrons with the spin-dependent chemical 
potential ]!„, defined by Eq. (fT2)) . Both Z" and S° are 
functions of the zero-Matsubara components ip^ and (pg of 
the bosonic fields, over which the integration in Eq. (|A1[) 
is carried out with 
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The functional integration over the remaining compo- 
nents of the bosonic fields results in the appearance of 
the phase correlation functions: 
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(A3) 



The functional averaging above is carried out with the 
weights exp[— 5'^'''], where S'^'^ are obtained from the ap- 
propriate bosonic action in Eq. ([5]) by subtracting the 
zeroth Matsubara components of Eq. (|A2p . 

The calculation of the the correlation functions of 
Eq. (|A3|) is carried out in Appendix [C] using the tools 
of stochastic analysis. The results are: 
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(A4) 



^ j^A^.n'-^a' =ns^^~,g2,.iJV-g 



At this point we introduce the winding numbers, as dis- 
cussed in section lTVl of the main text. Following the tran- 
sition (^Q = ujm + ^, and utilizing the identities Z.°(/i — 



iujm) = Z°{fj.) and S"{T,fi - iw,„) 
end up with 



-iuJrn'T GO 



g°(r,/x), we 






OO TT 



d^o/d^Se-^-^+-~°('"»-^^Z''(/I)S°, 



The grand partition function Z(/i) above is represented 
by the same double-integral and sum with S*' replaced 
by 1. The exponential factors involving N and n^ aris- 
ing from the phase correlation functions and the non- 
interacting GF cancel each other out exactly. This is 
hardly surprising as they are completely arbitrary. 

The summation over the winding numbers above can 
be performed using the Poisson formula 

^ /(27rfc) = - ^ / e""-/(^)dx. (A5) 



??l — — CXD 



This results in the expression for GF given by Eqs. (fT7|) 
and (J18p in the main text. 



Appendix B: Calculations in Auxiliary Canonical 
Ensemble 



We express Z°{^) in Eqs. ([T7)) and pS|) via the sum of 
the canonical partition functions for a system of n non- 
interacting electrons, Z,*^, using the standard relation 



To express the results of further integration in a conve- 
nient way, we also define the grand canonical and canon- 
ical partition functions with one level, Eq, excluded: 



a'y^a 



(B2) 



Then we substitute into Eq. ()18p the finite temperature 
GF of non-interacting fermions 

SlAr >0,^l)^ c-(^^°-^)- (1 - n„,.(M)) , (B3) 

where na_a{P') is the Fermi-Dirac occupation factor. We 
limit the calculation to Q{t > 0), since 5{t) = —9{t + j3). 
RecaUing that Z^iJi) = Z^(/It)Zj(/I|) we cast Eq. ^ 
into the form 



3n — 






Me'^^oe-(^°-^+^f)^ 



x\^ Qll3Km+n)-iipg{m+n)-a{m-n)ipl]<^ 11 \2^_ 



Carrying out the integration over (^q yields a Kronccker 
delta (57v,n+m- Performing the Gaussian integration over 
Lp% and defining M = m — n we find 



N 
M=-N 



fi^N+if}JiM+y,)^He^-^)r2^^^^^^Z, 



Substituting this into Eq. p7|) yields after straightfor- 
ward algebraic manipulations 
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^W N M=-N 
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where C7v,j\/ are defined in Eq. ([M|) . 

The canonical partition functions Z^r and Zjv( 
evaluated in Appendix |Dl resulting in 



(B4) 
to) are 



Z^. = e-^^^^^ 



ZnHo) = [1 - "^(mo)] Zw , (B5) 



where the Fermi factor for the a*^ level, na{^o) = 

[l -I- e^*^'^""'^"-'] , is taken with the auxiliary chemical 
potential j2o = A^A. 

Substituting Eq. jH]) into Eq. (1B4|, we find: 



N 
M=-N 



(B6) 



where /i and A'o are defined in Eq. (|25p and the auxil- 
iary partition function Z = F(0) is given by Eq. ([^0]) in 
the main text. The result is naturally spin- independent. 
Technically, the formal spin dependence vanished when 
calculating the integral 3n, Eq. (fT8|) . The GF for neg- 
ative T can be obtained from Eq. (jB6p using S(— t) = 
-S(/3-t). 

Now we find the full GF by summing over all sin- 
gle particle states Sa- This summation is carried out 
in the usual way by making the substitution ^^ Sa -^ 

OO 

A^^ / 9{ea)dea, i.e. effectively by averaging over disor- 



der by introducing the mean level spacing A. This leads 
to Eq. P^ in the main text. 

Finally, wc write the GF in the energy representation. 
Making the standard analytical continuation to the real 
time, T — > it, and Fourier transforming the GF to the 
energy domain we obtain the retarded GF used in the 
calculation of the TDoS as follows: 



(B7) 
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_e - 


- £« - 'Ca'.M + iO £-£„ -^Ar_i,A/-l +iO. 
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Appendix C: Phase Correlation Functions and 
Stochastic Analysis 



Now we write 11'^ (t), defined in Eq. (|A3p . in terms of 
the transition probability function (|C4p : 



Here we use stochastic analysis to calculate the phase 
correlation function 11'^ defined in Eq. (jA3[) . ff has been 
calculated in exactly the same manner. 

We note that the gauge equation (|10a|) can be viewed 
as a Langevin equation wherein the field ip'^ (r) plays the 
role of the stochastic force (noise), the distribution of 
which is governed by the action S^^ obtained from the 
appropriate bosonic action in Eq. ([5]) by subtracting the 
zeroth Matsubara components of Eq. (IA2[) . The noise 
correlation function is given by 



{^%tW{t')) = 2E, 



^(r-r')-l 



(CI) 



which follows from the expansion of 'p'^{t) in terms of 
Matsubara components: ^''(t) = X^m^^o "^m^"'""^- In- 
deed, the functional distribution of ip'^ is 



'D[ip''{T)]c- -^0 '^^i^"^^^^^^"^''^"'-''^} 
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(C2) 



which corresponds to {'fm'Pn) ~ '^EcP^^S^.-n, immedi- 
ately leading to Eq. (|Cip . 

It is convenient to represent the noise field as <^'^(t) = 
77(r) + iC with 77(r) a random function and ( a Gaussian 
random variable satisfying (^(r)} = {() = {'>i{t)0 = 0, 
ivirMr')) = 2EJ (r - r') and (C^) = 2E,/I3. As rjir) is 
standard white noise, we follow the standard procedure^^ 
to map the Langevin equation (jlOap to a Fokker-Planck 
(FP) equation: 
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Here CP5(0,t;0', r') is the conditional transition prob- 
ability function for a given (, formally defined by 

y^(6l, r; 9', r') = /s [e{T) -e]S [0{t') - 0'] \ where the 

(. . .)n means averaging over the white noise rjlr). The 
full transition probability function 3"^(0, r; 9', t') is given 
by the subsequent averaging over the quenched (i.e. r- 
independent) variable Q (as, e.g., in Ref . 1 171) : 

y^{9,T-e'y)^i5[e{r)-e]5[e{T')-e']) , (C4) 
i.e. ?^ ^ (y^^)^. 



Equation (|C3|) is a standard diffusion equation with a 
drift term. Its solution, with the natural boundary con- 
dition J"^(6', r; 9', t\C,) ~ 5{9 — 9'), is a decaying Gaussian: 



exp 
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= fd9d9'y%9,n;9',Tf)e-'^'^-'''K 



,i[e<^(r^)-r(T,)] 



Substituting here the solution (jC5p . we find the condi- 
tional (for a given C) phase correlation function as 

n^(r) = e-^-^e-^"e"''"^^ , (C6) 

where we defined t = |r/ — Ti|. Finally, the averaging 
over the quenched random variable C results in the first 
of Eqs. (|A4[) . The second one, for n'^(r), has obtained by 
applying, step by step, exactly the same procedure. 



Appendix D: Calculation of the Canonical Partition 
Function 



In this appendix we evaluate the canonical partition 
functions Zn and Z^^^a) defined in Eqs. (jBip and (|B2p . 
It follows from Eq. ^^ that 
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We calculate ZM{^a) (and thus Zm) in the saddle-point 
approximation: 



ZnI^c) «e-s=(^o) r '^c-k_[K(^l)]{^'-^l? ^ (D2) 
./_^ 27r 



where 

§o((p^) = -iiV^^ - In n (l + e-^'-'-"^°) . (D3) 

The saddle-point equation, §^ = 0, is convenient to write 
by replacing Y.a' fi^a') with A'^ /p°° de/(e) as in Ap- 
pendix (|B]). This gives, after calculating the integral, the 
following equation for finding (^[j: 



N- 



l + c'^^°+''^o /3A 



-— In 1 + e-''^" 



(D4) 



The Fermi- factor there, being of order 1, can be ne- 
glected, which means that the same saddle-point we 
would find in a calculation of Zj^: for large enough A^ 
the saddle-point is unaltered by the exclusion of a single 
state. Assuming also that N is so large that /3NA ^ 1, 
we find from Eq. (jD4p : 



(D5) 



-it^s = m^ = /3m- 
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In the same approximation Sa{(p'^Q) = l/(/3A), so that 
calculating the Gaussian integral in Eq. (|D2p gives 



Z 
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\PJZN 



(D6) 



while ZnI^o) differs only by the exclusion of the level a: 

^Ar(^a)-[l-n„(A2)]Zjv. (D7) 
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